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We study the CU5O4 cluster by exact diagonalization of a three-band Hubbard model and show 
that bound electron or hole pairs are obtained at appropriate fillings, and produce superconducting 
flux quantisation. The results extend earlier cluster studies and illustrate a canonical transformation 
approach to pairing that we have developed recently for the full plane. The quasiparticles that in 
the many-body problem behave like Cooper pairs are W=0 pairs, that is, two-hole eigenstates of the 
Hubbard Hamiltonian with vanishing on-site repulsion. The cluster allows W=0 pairs of d symmetry, 
due to a spin fluctuation, and s symmetry, due to a charge fluctuation. Flux quantisation is shown 
to be a manifestation of symmetry properties that hold for clusters of arbitrary size. 
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I. INTRODUCTION 

In many high-T c cuprates, one has superconductiv- 
ity at concentrations about 1.2 holes/Cu atom, that is, 
somewhat above the antiferromagnetic region at half fill- 
ing. In some cases, one finds superconductivity below 
half filling. Usually, one considers half filling as the vac- 
uum and speaks of electron superconductivity in such 
cases. Electron pairing is actually realized M, e.g., in 
the V structure of (Nd,Ce) 2 Cu0 4 . The T structure 
of this compound is different from the T structure of 
La2CuC>4, but is still characterized by CuC>2 planes [Q. 
However, increasing experimental evidence obtained in 
several cuprate superconductors suggests that the pairs 
exist above the critical temperature either in the form 
of superconducting fluctuations or preformed pairs. The 
latter aspect is apparent in the underdopcd (normal) re- 
gion in which a clear pseudogap essentially of the same 
magnitude as the superconducting gap is measured ||. 
The pairing state of these materials has d — wave sym- 
metry, probably mixed with s — wave |Q| . All these sig- 
natures put strict constraints to any microscopic model 
of the cuprates. Any theory of the paired state must 
predict d and s symmetries, and the pairing mechanism 
must be robust. It must survive well into the normal 
state, and operate in a wide range of concentrations far 
from optimum doping. 

In BCS theory, the first-order repulsion between like 
charges is overcome by the second order interaction 
with phonons. In high-T c superconductors the electron- 
phonon interaction is strong and phonons must be ex- 
pected to contribute in an important way to the pairing 
interaction, although their task looks harder because the 
repulsion integral U is large (several eV). However, the 
straightforward idea that the high-T c phenomena are just 
a rescaled version of BCS theory is not granted. The role 
of phonons may be important, but is different, and some 
other ingredient is essential. 

Rather than proposing a new attractive interaction, we 
wish to point out that the popular repulsive Three-Band 
Hubbard model already leads to pairing. The model is: 



H = Hn 



w 



(1) 



where the independent hole hamiltonian reads, in the site 
representation 



Cu 



(2) 



where n.n. stands for nearest neighbors ||. The on-site 
repulsion Hamiltonian will be denoted by 



(3) 



where Ui — Ud for a Cu site, Ui — U p for an Oxygen. 



As in previous work |15|, we use standard parameter val- 
ues, i.e., Ud=5.3,U p = 6, t=1.3, e<j=0, e p —i.b. The hole 



parameters Ud = 5.3 eV, U p — 6. eV differ somewhat 
from other literature estimates fltifl , and must depend 
on the compound and doping. For La2Cu04, U p —4 eV 
and £7^=10.5 eV have been recommended ||. None of 
our results depends qualitatively on the precise value of 
the model parameters, since ours is basically a symmetry 
argument. The electronic properties of this model are un- 
der intense investigation by several approximations based 
on perturbation theory and the Bethe-Salpeter equation. 
The FLEX approximation is a generalized RPA |7) and 
leads to pairing and superconductivity in the three-band 
Hubbard model |g. The excitation spectra of the 2D 
Hubbard model have also been studied by a related self- 
consistent and conserving T-matrix approximation by 
Dahm and Tewordt ||; we mention incidentally that 
recently diagrammatic methods have been successfully 
applied to the photoelectron spectra of the Cuprates in 
other contexts too, like the spin-fermion model |[(^. A 
perturbative expansion around the strong coupling limit, 
in powers of the kinetic energy, requires a nonstandard 
cumulant expansion, but is feasible, as shown quite re- 
cently by Citro and Marinaro [|ll| for the p — d model 
which is the present model with U p = 0. In this way, they 
have shown that normal state properties like the specific 
heat as a function of doping can be well understood |L2j ; 
they also derived the effective pairing interaction in the 
same approximation p3| and studied the doping depen- 
dence of the superconducting transition temperature jl4| . 

Our starting point is the observation that, due to the 
planar symmetry, there is actually no repulsion bar- 
rier to overcome. In a series of papers Jl5| , p7[p^ ] we have 
introduced the two-hole singlet eigenstates of the Hamil- 
tonian with zero Coulomb on-site repulsion (the so called 
W=0 pairs). They arise in the full plane, and also in the 
clusters that possess the same full symmetry around 
a central Cu as the full plane. In the full plane, this 
situation is always realised, because W=0 pairs can al- 
ways be obtained from holes at the Fermi level; in clus- 
ters, on the other hand, the hole number (relative to 
the true hole vacuum) must be such that two holes par- 
tially fill a degenerate state. In the many-body problem, 
two holes at the Fermi level in a W=0 pair state do not 
interact directly; however the pair is dressed by the in- 
teraction with the background particles. By exact diag- 
onalisation of cluster Hamiltonians with up to 21 atoms 
and 4 holes, we demonstrated [p| that the dressed W=0 
pair is a bound Cooper pair, and quantizes the magnetic 
flux like superconductors do. Any strong distortion of 
the cluster symmetry breaks the pairing and restores the 
normal repulsion JL7J. We also considered first-neighbor 
0-0 hopping and off-site interactions p7j . Remarkably, 
the off-site repulsive interactions, when included, tend to 
enhance the effect somewhat, so we devote the present 
study to the on-site interaction effects. The binding en- 
ergy of the pairs in these clusters is of the order of tens of 
meV, which is not comparable to any of the U and t in- 
put parameters. The reason is that the interaction, which 
vanishes identically for the bare W=0 pairs, remains dy- 
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namically small for the dressed quasiparticles. Indeed, by 
a diagrammatic analysis we demonstrated that low-order 
perturbation theory is a good approximation to the exact 
diagonalisation results and allows to understand that the 
attraction in the d channel is due to virtual spin-flip ex- 
citations. This suggests that a weak coupling theory may 
be useful to study the pairing force, despite the fact that 
U is not small compared to t. It is an obvious limitation 
of the cluster approach that W=0 pairs are possible at 
discrete values of the hole concentration. Our previous 
cluster calculations suggest that the mechanism operates 
in a much broader range of hole concentrations than is re- 
alised in actual Cuprates, from very highly overdoped (as 
in CuOi) to very low (as in Cu^Oiq). The diagrammatic 
analysis further demonstrates that the effective interac- 
tion is the result of a partial cancellation of positive and 
negative contributions, so it is not necessarily attractive 
in all cases; the general signature of W=0 pairs is that 
the absolute value of the interaction is much smaller than 
in the other cases. 

The mechanism we are considering is only a part of 
the story, but it seems to be a most peculiar part, being 
related to nothing but the Ci V symmetry. For similar 
reasons here we wish to make abstraction from phonon 
effects to see how far the idealised description can account 
for reality by itself. We believe that a mechanism which 
predictably gets attraction out of repulsion is by itself of 
theoretical interest. 

Next, we have generalised the theory of pairing to the 
full plane ppf . In short, one finds W=0 pairs at the Fermi 
level for any concentration and this leads to a Cooper-like 
instability of the Fermi liquid. Pairing prevails for a range 
of concentrations above half filling, in agreement with the 
results jjij of the Renormalization Group technique. We 
have shown that the full configuration interaction calcu- 
lation can be performed recursively. At each step, one 
decouples a class of virtual excitations while renormal- 
ising the matrix elements of Hq and W . At the end, 
one obtains an exact, analytical canonical transforma- 
tion producing an effective Hamiltonian for the dressed 
pair. In order to get actual numbers, however, we had 
to neglect the renormalisations in the final formula; this 
approximation is fully justified at weak coupling. 

In the present paper we extend the analysis of Ref. |lg| ] 
by diagonalizing the CU5O4 cluster with increasing num- 
ber rih of holes. We wish to demonstrate that the power 
of the symmetry driven mechanism is such that attractive 
interactions arise even in small clusters with > 4, de- 
spite the high hole concentration. One can proceed from 
the true hole vacuum and insert holes until the last two 
form a W=0 pair; if the interactions produce a bound 
state we conventionally speak of hole pairing. Alterna- 
tively, one can proceed from the true electron vacuum 
and insert electrons until the last two form a W=0 pair; 
if the interactions produce a bound state we convention- 
ally speak of electron pairing. These two expressions sim- 
ply mean that we get pairing by adding two holes (as in 
La2Cu04) or two electrons (as in (Nd,Ce)2Cu04), re- 



spectively |22j. The physical point here is that electron 
pairs and hole pairs are related by a charge conjugation 
symmetry and the very same basic mechanism or dia- 
gram is operating in both cases. We find new instances 
of electron and hole pairing, again with a binding energy 
of a few tenths of meV in the physical parameter space. 
We also find a case when a W=0 pair leads to a weak re- 
pulsion. Further, we demonstrate how the two different 
symmetries (Ai and B2) of W=0 singlet pairs allowed by 
the cluster are both necessary to produce the supercon- 
ducting flux quantisation phenomenon in this cluster. 

One reason for considering clusters again, after much 
excellent work from several authors p3| and our own pre- 
vious work on clusters and on the plane, is that in this 
way we can test our canonical transformation approach 
against the results of exact diagonalisation. Another rea- 
son is that we wish to explore the relation of the flux 
quantisation phenomenon to the symmetry Group in the 
presence of the vector potential, which breaks the transla- 
tional symmetry. Our main questions are: can pairing be 
reliably predicted by studying the behavior of the system 
at weak coupling? Is the superconducting flux quantisa- 
tion property exclusive of small clusters, or is it a general 
consequence of symmetry? 

II. ONE-BODY ENERGY LEVELS OF THE 
SYMMETRIC 9-SITE CLUSTER 

In the hole picture, the one-body energy levels of 
CU5O4 are those displayed in Table |. Here, the hole 
vacuum is a state with no holes at all. 

Two levels are triply degenerate, comprising twice de- 
generate states of e(x,y) symmetry and states belonging 
to 61; this accidental degeneracy is due to the fact that 
in this small cluster any permutation of the four Cu — O 
units bound to the central Cu is a symmetry; therefore, 
the full symmetry Group of the cluster is S4, which has 
Ci v as a subgroup, and admits degeneracy 3. Since this 
property does not extend to the plane, we continue using 
the irreducible representations (IRREPS) of C± v anyhow. 

In the electron picture, the levels are met in reverse 
order, but the sequence of symmetry labels remains the 
same. Thus, one notices that there is an approximate 
electron-hole symmetry, or charge conjugation symmetry, 
in this model. 



III. W=0 PAIRS 

Both in the full plane and in clusters, the W=0 pairs 
are due to the symmetry, but there are some differences 
between the two cases, that we wish to stress in this 
Section. Let us first review the theory for the plane [^0| . 
Omitting the band indices, we shall mean 

\d[k]) = \\k + ,-k.\\=cl + ci kt _\vac> (4) 
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to be a two-hole determinant al state derived from the 
Bloch eigenfunctions (\vac > is the true hole vacuum). 

The point symmetry Group of the Cu-0 plane is C\ v . 
We introduce the determinants Rd[k] = d[Rk], R £ C^ v , 
and the projected states 
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^[ fc ] = 7f E x (r,) (R)\Rd[k}) 



(•5) 



R£C iv 



where is the character of the operation R in 

the Irrep rj. In the non-degenerate Irreps, the opera- 
tions that produce opposite Rk have the same character, 
and the corresponding projections lead to singlets. Let 
Ri,i = 1, ..8 denote the operations of C± v and k, kl any 
two points in the Brillouin Zone (BZ). Consider any two- 
body operator O, which is symmetric (RjORi = O), and 
the matrix with elements Oij —< d[k]\Rl(JRj\d[k'] >, 
where k and k' may be taken to be in the same or in 
different bands. This matrix is diagonal on the basis of 
symmetry projected states, with eigenvalues 

0( Vl k 1 k') = Y / X {7 HR)0 R (k,k') (6) 



R (k,k') = {d[k]\6\Rd[k']). (7) 



where 



Thus, omitting the k, k! arguments, we get in particular 

O (%) = O e + C2 + Ci + O c s 

-O a , - O ay - O a . - O a , (8) 



O ( 1 B 2 ) =0 E + C2 - 

-o a . - 0„ 



o Ci - 
+ o„, 



-o„, 



(9) 



If O is identified with W, since W E = Wc 2 = = W ay 
and W Ci = Wc 43 = W a[ = W^, one finds W ( 1 A 2 ) = 
W ( 1 B 2 ) = . These are W=0 pairs, like those studied 
previously fisf in clusters. In the full plane, however, 
W=0 pairs are obtained from holes at the Fermi level for 
any filling. 

Small clusters like CU5O4 allow a nice illustration of 
the theory because they also allow W=0 2-body solu- 
tions. This property is a consequence of their full C± v 
symmetry around the central Cu. However, there are no 
Bloch states in a finite cluster with open boundary con- 
ditions, and the W=0 singlet pairs come out differently 
First, we may consider the orbitals of (x,y) symmetry in 
Table and form 2-hole determinants 



d[x,y] = ||ar+,y_|| = cl + cl_\vac > 
d[y, x] = \\y + ,x-\\ = 4 + 4_ \vac >; 



they are eigenstates of H and have the W=0 property, 
since the amplitude of double occupation of any site is 0. 



Unlike the case of the full plane, no projection like that 
performed in Equation (Q) is necessary here to get the 
property. In CU5O4, there are two sets of (x,y) states, 
so the x and y of the above equation may belong to the 
same or to different sets. If the x and y states are taken 
from the same set, the singlet 



^B 2 ) = 



d[x, y] + d[y, x] 
V2 



(11) 



is an eigenstate of the kinetic energy and of W, belongs 
to the W=0 eigenvalue and to 1 B 2 . W=0 pairs of this 
symmetry and of l A 2 exist in the full plane as well. If 
the x and y states are taken from the different sets, we 
denote one of the sets by a prime and consider two-hole 
determinants like d[x, y']; these are eigenstates of H and 
are W=0 pairs, however they do not belong to any of the 
IRREPS of C± v . We can form singlet combinations with 
the W=0 property both in the 1 B 2 and l A 2 symmetry, 
namely, 



d[x, y'} + d[y, x'} + d[x', y] + d[y\ x] 



which belongs to 1 B 2 and 



d[x, y'} - d[y, x'} + d[x', y] - d[y', x] 



(12) 



(13) 



which belongs to 1 A 2 . In addition, there are also 1 Ai 
W=0 pairs, using the degenerate x, y and b = b\ orbitals. 
The two-hole states 



\X+X- 



\y+y- 



V2 



\x 2 +y 2 > 



and 



16+6-1 



\bb > 



(14) 



(15) 



are a basis of degenerate eigenstates of Ho having x Ai 
symmetry. Diagonalising W in this basis we get two-hole 
eigenstates of H . The 2x2 matrix of W is: 



<bb\W\bb> < bb\W\x 2 + y 2 > 

< bb\W\x 2 + y 2 > < x 2 + y 2 \W\x 2 + y 2 > 



Up 

4 

U P 

2V2 



U P 

2J2 



The lowest eigenvalue is and the W=0 pair is 



\\x l + <f 



^A 1 ) = -\l-\bb) 



(16) 



(17) 



This type of W=0 pairs does not exist in the full 
(10) plane. The upper eigenvalue is —f- and the eigenfunction 



\\bb) + ^Jl\x 2 +y 2 > is strongly affected by the on-site 
repulsion. 
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Now consider the C115O4 cluster in the non-interacting 
limit with 2 holes, which sit the lowest level of a\ sym- 
metry that we denote a for short. Let this be the new 
vacuum state |0 >. Adding 2 holes, we partially fill the 
next degenerate levels, which can give rise to W=0 pairs. 
Starting with the above defined two-hole d determinants, 
we can form 4-hole determinantal states, denoted by a 
capital D, like for example 

D[x,y] = cl + cl_\0> 

= \\x+a+,y-a>-\\i ( 18 ) 

in keeping with the notation of Ref . |2(J , below we shall 
denote such configurations with the background a or- 
bital occupied for both spins as m states. We shall also 
need the other 4-hole determinantal states, namely, the 
a states, in which one of the background a spin-orbitals 
is not occupied by holes, and the (3 states in which both 
background a spin-orbitals are missing. For the C115O4 
cluster in the non-interacting limit a similar situation is 
realised if the new vacuum state is taken with 10 holes, 
filling the lowest levels according to the aufbau principle. 
Adding 2 holes, again we partially fill the next degener- 
ate levels, which can give rise to W=0 pairs. A similar 
definition of m, a and (3 is possible, and more excited 
configurations also exist. 

The matrix elements of W in this model have no ex- 
change terms, since only holes of opposite spin can inter- 
act. The diagonal elements W m ^ m can be expressed in 
terms of orbitals p, q, r, s by the two-hole integrals 

W(p,q,r,s) =Y,U lP *(i)q*(i)r(i)s(i). (19) 

i 

For example, the m state of Equation(|l^) yields 

W m ,m = W(x, a, x, a) + W(y, a, y, a) + W(a, a, a, a). 

(20) 

We note that this is the Hartree-Fock interaction. The 
last term of the expression refers to the interaction be- 
tween the background particles in the a spin-orbitals and 
is the same for all the m states, and the rest brings out 
single-particle corrections to the energy of the orbitals 
and could be readsorbed in the definition of Hq. The im- 
portant point is that no term contains both x and y; no 
direct interaction between the two added particles exists, 
because of the W=0 property. 

The matrix element of the two-body operator W be- 
tween determinants which differ by two spin-orbitals are 
given by the well-known rule 

(\\k + ,k-,ui...Un\\W\\k' + ,k'_,ui...u n \\) 

= (\\k + ,k_\\W\\k' + ,k'_\\) (21) 

where k is different from k' while u\ . . . u n is a sequence 
of occupied spin-orbitals. Using Equation (0), one finds 
that in the full plane, the matrix elements between dif- 
ferent m states W m , m > vanish. In the cluster, this is not 



true. Let x denote the state of x symmetry taken from 
the lower e degenerate level, and y' denote the state of y 
symmetry taken from the upper e degenerate level. Then, 
the m states involving these orbitals are coupled by W 
to those involving x and y, and to those involving x' and 
y' . Such matrix elements are forbidden in the full plane 
by momentum conservation, but exist in finite systems 
with open boundary conditions, having no translational 
invariance. They couple m states belonging to different 
eigenvalues of H . 

IV. THE EFFECTIVE INTERACTION 

We need a rigorous definition of the effective interac- 
tion between two holes in many-body systems, and this 
requires a careful analysis. Actually, we shall use two 
alternative definitions, one of which is suitable for nu- 
merical exact diagonalisation work, while the other one 
is much more microscopic and analytical. Therefore, we 
have to show that these two definitions essentially agree 
and lead to the same physical conclusions. That will 
emerge from the analytical treatment of the present Sec- 
tion and from the numerical results of Section 6. 



A. First definition: A 

When working by exact diagonalisation, we consider 
a cluster with holes; its interacting ground state en- 
ergy Eh(jih), obtained with the Hamiltonian of Equations 
ffljplpj]), is referenced to the hole vacuum for any n^. In 
terms of these eigenvalues we define, following references 
(f§, ©), 

A h (n h ) = E h (n h ) + E h (n h - 2) - 2E h (n h - 1). 

(22) 

^■h(n-h) is one definition of the pairing energy. This def- 
inition is simple, but requires computing the eigenval- 
ues with great accuracy, and has several drawbacks. It 
says nothing about the dynamics which leads to pairing. 
Moreover, generally a negative A does not unambigu- 
ously imply pairing, and further problems arise |2l| since 
the above definition depends on the comparison of sys- 
tems with different rih- 

However, the application of Equation ( |22] ) is safe in 
the specific case when the last two holes are in a W=0 
state; in Ref. [j^), we have shown that in this case A 
really coincides with the ground state expectation value 
of the effective interaction, at least at weak coupling; if 
the interaction is attractive and produces a bound state, 
A/ l (n? l ) is negative and |Aft,(n/j)| is the binding energy. 
These results were obtained by analysing exact diagonal- 
isation results for clusters with nh = 4 by lowest-order 
perturbation theory. In the present paper, we wish to ex- 
tend those results to larger by exact diagonalisation 
and a more powerful analytical method. 
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B. Second definition: W, 



iff 



The alternative definition is intrinsic to the rih holes 
system and much more transparent. We achieve it by 
a canonical trasformation that determines the effective 
two-body Hamiltonian H from the many-body H of 
Equation (]l|). We set up the Schrodinger equation for 
the ground state of the cluster with holes, namely 



H\*o) = £o|*o) 



(23) 



Here, E = E^nh)- We take the ground state configu- 
ration of the noninteracting rih ~ 2 system as our vac- 
uum state (the non-interacting Fermi sphere). The exact 
l^o > can be expanded in terms of excitations over the 
vacuum: 

|*o >= E° m l TO > +Y^b a \a > + E c /3^ > +.... 

m a ft (24) 

here m runs over pair states, a over 4-body states (2 
holes and 1 e-h pair), (3 over 6-body ones (2 holes and 
2 e-h pairs). In C115O4 with 4 holes, the vacuum is the 
a\ configuration and the expansion terminates with the 
f3 states; it terminates anyhow in any finite system, after 
a finite number of terms, so there are no convergence 
problems. Next, we consider the effects of the operators 
on the terms of |*o >• We write: 



H Q \m >= E m \m >, H \a >= E a \a >, . 

and since W can create or destroy up to 2 e-h pairs, 
W\m >= ^2W m ', m \m! > +El a > W a , m 

m' a 





(25) 



(26) 



For clarity let us first write the equations that include 
explicitly up to 6-body states; then we have 

W\a >= ^2\m > W m , a + Y)c/ > W a ,, a 

m a' 

+Y)P > Wp, a (27) 


where scattering between 4-body states is allowed by the 
second term, and 

W\P >= E l m '> W ™',f> + E I") W «,0 



W) w >, 

0> 



(28) 



The Schrodinger equation (|23J) yields equations for the 
coefficients a,b and c 



(E rn — E ) 

+E a ™' W m , m > + ^2b a W m , a + E^^W = (29) 

m' a ft 

(E a — Eo)b a 

+ 1 £ / a m >W a , m > + Yy^' W ^' + ^2 c /3W a>0 = (30) 

m' a' ft 

(Eft — E )cft 

+E°W Wft. m > + J2 b *'Wft, a > + J2 c 0' W 0-0' = ( 31 ) 

m' a' ft' 

where Eq is the interacting ground state energy. In prin- 
ciple, the Wft^ft term can be eliminated by taking linear 
combinations of the complete set of j3 states. The com- 
plete set of (3 states can be chosen in such a way that 



(H + W)ft,p, =E' S(p 1 p') 



(32) 



with this choice, the Wp^p terms are removed, while 
E'a replaces the noninteracting eigenvalue Eft. In other 
terms, we get a self-energy correction to Ep and a mix- 
ing of the vertices, without altering the structure of the 
equations. Then, we may rewrite Equation ( |3l| ) in the 
simpler form 



(E'p - E ) eft + ^a m >Wft, m > + ^b a >Wi3,a> = 0. 

m' a 1 (33) 

Now, we exactly decouple the 6-body states by solving 
the equation (|33] ) for cp and substituting into ( p9l^p| ), 
getting: 



(E m — E Q ) a m + y^a,, 

+E 6 « 

a 

(E a - E ) b a + Y, 

m ' 



w. 



"' " ! r 7 ■ TP TP' 







w„ 



E 



W^ftWft^ 
'0 - "p 



En — E Q 



(34) 



W a ,m' + E 



W a ,ftWft, 





W ^ + ^E -E> 



Eq - Ep 

w atf3 Wft ta 1 











(35) 



Introducing renormalised interactions W , we may 
rewrite these equations in the form 



(E m - E ) a m + E«m'W4, m , + E^^C* = 



(36) 
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(E a — Eq) b a 



]b a >W^ a , =0 



(37) 



If in Equations (|2^,|3(]) we drop the terms involving 
the f3 states, they reduce to the same form as Equations 
( [36],[57]), except that in the latter equations some quanti- 
ties are renormalised. In other terms, the role of 6-body 
states is just to renormalize the interaction in the equa- 
tions for the 2-body and 4-body ones, and for the rest 
they may be forgotten about. If E is outside the con- 
tinuum of excitations, as we shall show below, the cor- 
rections are finite, and experience with clusters suggests 
that they are small. Had we included 8-body excitations, 
we could have eliminated them by solving the system for 
their coefficients and substituting, thus reducing to the 
above problem with further renormalizations. This is a 
recursion method to perform the full canonical transfor- 
mation; it applies to all the higher order interactions, 
and we can recast our problem as if only 2— and 4-body 
states existed. 

Again, the W a , term can be eliminated from Equa- 
tion (p7|) by taking linear combinations of the a states. 
This is achieved by choosing the complete set of a states 
in such a way that 



{H + W') c 



E'J(a,a'). 



(38) 



With this choice, the W' a a , terms are removed, while 
E' a replaces the noninteracting eigenvalue E a . In other 
terms, we get a self-energy correction to E a and a mix- 
ing of the vertices, without altering the structure of the 
equations. Now Equation (B^) becomes 



(E' a -E )b a +J2arn'W^ m , =0 



(39) 



Solving Equation ( p9| ) for b a and substituting in Equation 
( |36| ) we exactly decouple the 4-body states as well. The 
eigenvalue problem is now 



(E - E m ) a m = 2J a m ' (m\S[E ]\m') 



(40) 



where 



\S[E ]\m')=W^ m ,+J2 



< m\W'\a >< a\W'\m' > 



Eq-EL 



(41) 



We introduce the diagonal elements of the a summation: 
p _^2 <m\W'\a><a\W'\m> ^ ^ 



En-E' 



then, Equation (|4(]) becomes 



Eaa m — [E„ 



where for m ^ m' 



W'+E 



m\W eff \m') 



(43) 



/ ,Txr , a rxr/ ^{m\W'\a){a\W'\m') 
(m \W eff \m') = W' mtm , + }2— — 



Eo~E' a 



(44) 



The W m m , term does not arise in Ref. |2(| because in 
the full plane it vanishes by momentum conservation. 

Equations ( p^|]44| ) determine the amplitudes a m of the 
m states in the n^-hole state and the ground state eigen- 
value Eq relative to the hole vacuum. Their solution, 
inserted in Equation ( |39|) yields the coefficients b a and 
we could proceed with the full calculation of ^>q; this ap- 
pears to be hard for a large system. However, our task 
here is to find the effective two-body Hamiltonianian; this 
is much less expensive. 

Indeed, Equation ( (43| ) is of the form of a Schrodinger 
equation with eigenvalue Eq for pairs with effective in- 
teraction We//. Then we may interpret a m as the wave 
function of the dressed pair, which is acted upon by 
an effective Hamiltonian H. The change from the full 
many-body H to H is a canonical transformation which 
holds to all orders. W e // is the effective interaction be- 
tween dressed holes, while F is a forward scattering op- 
erator, which accounts for the self-energy corrections of 
the one-body propagators: it is evident from ( (43|) that 
it just redefines E' m . Also in Cooper theory |29| one 
meets electron-phonon self-energy terms, which do not 
contribute to the effective interaction. The basic spin- 
flip diagram responsible for We// had been identified 
before pq |. Any other pairing mechanism not consid- 
ered here, like off-site interactions, inter-planar coupling 
and phonons, can be included as an extra contribution 
to W m , m which just adds to W e //. 

This way of looking at Equation (|4^) is perfectly con- 
sistent, despite the presence of the many-body eigenvalue 
Eq, because we are not compelled to reference the energy 
eigenvalues to the hole vacuum. We note that if we shift 
Hq by an arbitrary constant AE in Equation (|25|), by 
setting 



H'=H a - AE 



(45) 



the same shift applies to the eigenvalues E m , E a , Ep and 
so on, and also to the renormalised quantities like E' a , E' g . 
Therefore, the effective interaction W e / / of Equation ( |44| ) 
and the F matrix elements are unaffected by the shift. 
Thus we can reference Eq to a new energy origin by shift- 
ing the diagonal terms in Equation(|43"|) without changing 
the off-diagonal terms. Since we wish to regard Equation 
( [43| ) as a Cooper- like equation for the pair, it is natural 
to set AE equal to the interacting ground state energy 
eigenvalue for the nt — 2 hole system, relative to the hole 
vacuum, 



AE = Eh(rih — 2) 



(46) 
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This quantity is obtained by diagonalising the cluster 
Hamiltonian with nh — 2 holes. In Ref . J2(]] , dealing with 
the infinite plane, this was our choice. The energy of two 
independent holes, relative to the rih — 2 background, is 
2Ep, where Ep is the Fermi energy; when the effective 
interaction is accounted for, the energy of the two bound 
holes is 2Ef + A, where |A| is the binding energy. 

Up to this point, the treatment is exact. However, we 
can make an easy use of Equation (|43] ) if we can neglect 
the renormalisations in Equation (|44|), setting W' — > W 
and E' a — ► E a , which is fully justified in the weak cou- 
pling case. This is the approximation that we proposed 
in Ref ( |2(J) and that we want to test in the present 
paper. In fact, if we are primarily interested in the sym- 
metry of the ground state, and in the presence or absence 
of pairing, we can get these results without a large com- 
putational effort. We exemplify the procedure for CU5O4 
in the rih — 4 case. Two degenerate to states are lowest 
in the non-interacting limit, namely, the configuration 
to = D[x,y] of Equation ( [l8| ) and to' = D[y, x], where 
the x, y orbitals belong to the lower e level; the ijji} A{) 
state also is degenerate with to, to', but by symmetry 
W cannot mix it to them. As already noted, m and to' 
do not interact through the W m m > term. To calculate 
W e ff, we rewrite W (Equation ([|)) in the orbital repre- 
sentation, with 

orb 

4 = £<i|">4 (47) 

V 

where v runs over all the orbitals, obtaining 



W 



orb 

£ W(fi,p,p,a)cl t+ cl_i 



(48) 



The pair (p + , ct_ ) which is annihilated may correspond 
to (y+,X-), (y+,a_), (a+,X-), (a+,a_). The first choice 
gives nothing since it corresponds to a W=0 pair; the last 
choice yields a /3 state. To lowest order, only the a states 
contribute, involving the excitation of either the a + or 
the a_ hole. Many of the W matrix elements vanish 
by symmetry; we are going to neglect those connecting 
to excited x',y' orbitals, which occur at higher energies, 
because we are considering weak coupling. 

Considering the contribution of (y+, a~), one finds that 
the only a states coupled to D[y, x] by W are those of the 
form |/xaj/a;| = |/i+a+j/_ af_|, in which the hole in a_ is 
promoted to y_ while y+ is scattered into Therefore, 



Ecx — S a ~\~ S i 



where s y — e x . One finds 



and 



{\nayx\ W \yaxa\) = -W(y, a, /it, y) 



(\xaya\ W \/j,ayx\) — W(x, a, 11, x). 



(49) 
(50) 
(51) 



Therefore, taking into account that each of the two back- 
ground a holes can be promoted and this brings a factor 
of 2, using we obtain 



(m\W eff \m') 



^W(y,a,fi,y)W(x,a,n,x) 



E a 



(e a + £ p + 2e x ) 



(52) 



Since W(y, a, x, y) — W(x, a, y, x) — 0, the empty (of 
holes) states \i belonging to the e representation yield 0. 
The empty orbitals that contribute are those of the a\ 
symmetry, that will be denoted by a', and those of bi 
symmetry that we shall write b. 

The a' orbitals contribute to the repulsion, and the 
b orbitals to the attraction. Indeed, W(x, a, a', x) = 
W(y, a, a', y), and the contribution of the states of a,i 



symmetry is — 2^ a , 



since E a > Eq 



E -(e a +e a ,+2ex 

this is positive. On the other hand, W(x,a,b,x) — 
—W(y,a,b,y), since the orbitals of of b\ symmetry 
change sign for a ? rotation. Therefore the contribu- 
tion of the b states is attractive. This is an example of 
the interference of opposite contributions to W e ff, that 
we emphasised in Ref pfl| . 

Therefore, at this stage, a self-consistent treatment of 
Equation pg| ) must be sought, because W e ff depends 
on the eigenvalue Eq. A straightforward recursion ap- 
proach leads to a continued fraction solution, which has 
contributions from all orders of perturbation theory. 



C. Equivalence of the two definitions at weak 
coupling 

For weak coupling, however, a cruder but simpler ap- 

poximation is justified: one calculates We// neglecting 

all self-energy corrections, in such a way that Eq in 

(|52|) reduces to 2e x + 2e a \ this is the lowest (second- 
— (2) 

order) approximation W^A. In the same, lowest-order 
approximation, the shift in Equation (^), which recasts 
Equation ( pig] ) as a two-body problem, in the Cooper- 
like form, is AE — 2e a ; further, one considers only the 
mixing of the degenerate configurations to = D[x,y] and 
to' = D[y 7 x] . In the resulting 2x2 problem, the diagonal 

(2) 

entries are identical, and the J * interaction produces 

the off-diagonal elements, with the result that the singlet 

(2) 

is stabilised by and the triplet is destabilised by 

the same amount. Therefore §6|, A(4)( 2 > = wf) f . One 
obtains for A(4)( 2 ) the following second-order expression: 



-2 



E 



W(a, b, x, x) 2 
(e& - e a ) 



E 



A(4)< 2 > = 
W(a, a', x, x) 2 
{(-a' - e a ) 



(53) 



where the sums run only over the one-body states of a 
and b symmetry. This agrees with the result that we ob- 
tained earlier mm from a diagramatic analysis of Equa- 
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tion Q22|), Thus, the two definitions of the effective in- 
teraction lead to the same result, at least in the weak 
coupling limit. 

We stress again that the sign of A(4)( 2 ) is determined 
by the relative weights of the virtual excitations to the 
empty states of different point symmetry. In general, 
W e ff can produce attraction or repulsion, depending on 
the hole concentration. This crude approximation will 
turn out to be sufficient for qualitative purposes, i.e, to 
predict when pairing occours. 



V. W=0 PAIRS AND CHARGE CONJUGATION 

Consider the cluster with rih holes. As shown above, 
the interesting situation arises when rih is such that, fill- 
ing the levels according to the aufbau principle, the last 
two holes go to a degenerate level. Accordingly, we ex- 
pect that A^(n^) measures the effective interaction be- 
tween the holes of the W=0 pair. In Ref. fL8[| , we have 
shown that this is the case at weak coupling. If the inter- 
action is attractive and produces a bound state, |Ah(nh)| 
is the binding energy. This situation can be realised with 
rih — 4 in highly symmetric Cu-0 clusters containing up 
to 21 atoms The last 2 holes then go to the lowest 
level of e symmetry. 

According to Table I, the CU5O4 cluster has an upper 
e level, which is reached with 12 holes, so we are inter- 
ested in A/j(4) and A/,,(12). Moreover, we can exploit 
the approximate electron-hole symmetry of the problem 
to obtain two more interesting cases. The approximate 
symmetry consists in the fact that the same sequence of 
symmetry labels is obtained by reading Table [l| from up 
to down and in reverse order. The reverse order corre- 
sponds to adopting the electron picture and starting from 
the electron vacuum. Going to the electron picture, the 
three-band Hubbard Hamiltonian becomes: 

H = ^(2e, + Ui) - ^(si + Ui)a\ a a ta 

i ia 
- ^ t ij a ia a j<r + ^ UiTH+nj-, (54) 

<Z,J>(7 i 

W=0 electron pairs are obtained for n e =4 and 12 elec- 
trons. Letting now E e {n e ) denote the ground state en- 
ergy of the cluster with n e electrons, the effective inter- 
action between the two electrons in the pair is measured 
by 



A e (n e ) = E e (n e ) + E e (n e - 2) - 2E e (n e - 1). 



(55) 

Since the dimensionality of the one-body basis is 18, 

A e (4) = £ e (4)+£ e (2)-2£ e (3) = 
E h (U) + E h (l6) - 2E h (15) = A h (16) (56) 

and, similarly, A e (12) = A h (8). 



We recall that we speak of electron pairs when two 
added electrons partially occupy a degenerate state and 
of hole pairs when the same situation is reached by adding 
two holes; however the final situation is exactly the same. 
For example, consider the W=0 pair state of Eq. jll]). 
One readily verifies that in a canonical transformation 
from holes to electrons, putting a\ a —Ci„^ the two-hole 
state becomes a two electron state of the same form. 
Therefore the two-body W=0 state is invariant under 
charge conjugation, and if holes are paired, electrons are 
also paired. In order to avoid switching all the time be- 
tween the two equivalent pictures, below we discuss ev- 
erything in terms of holes. Summarizing the results of 
the present Section, we can test the effective interaction 
in CU5O4 by calculating A^(n/j) with — 4,8,12 and 
16. 



VI. NUMERICAL RESULTS AND DISCUSSION 

By an enhanced Lanczos routine, we computed the 
ground state of the CU5O4 cluster with even numbers 
ri/i of holes and vanishing z component of the total spin; 
the parameter values are specified in the Introduction. 
The maximum size of the matrices (15,876) occurs for 
rih = 8; the =12 and 16 cases are handled by trans- 
forming to the electron picture. The results for A/ 1 (n/ 1 ) 
are summarised in Table O. One sees that A;,, for = 
4,8,12 and 16 is much smaller in absolute value than 
for other fillings, as expected. This confirms that W=0 
pairs are involved. In particular, for = 4,8 and 16 
A/ 1 (n? l ) < and pairing occurs, while for nh —12 a small 
repulsion prevails. When pairing is obtained, this means 
that the renormalisation of the parameters inherent in 
the canonical transformation does not have important 
consequences. To see if the behaviour at =12 is an 
exception, we have repeated the calculations with scaled 
U values. Using U p = .06 eV and Ua — -053 eV , which 
are < < t and allow applying perturbation theory, we still 
get a positive result, namely Ah(12)=0.0034 meV. Thus, 
even second-order perturbation theory would suffice to 
predict A?,, (12) > in this case. 

For rih = 4 we have an analytic second-order result 
(Equation |53| ) and we can check its degree of validity by 
comparing with the exact diagonalisation values of A. 
Making use of the above standard values of the param- 



eters we calculate the relative error S 
turns out that S < .07 up to U/t wl 



A-A< 2 > 



It 



A+A( 2 ) 

Thus, already 
the second-order approximation is remarkably accurate 
in estimating the effective interaction. We conclude that 
our treatment based on the unrenormalised formula of 
Equation ( |i4| ) correctly predicts the presence or absence 
of pairing, depending of the hole concentration, and even 
a simple, second-order approximation to it has a semi- 
quantitative accuracy when compared with exact results. 

The indications that we may draw from this Section 
are: i) our mechanism is operating for a wide range of 
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hole concentrations and produces a much reduced inter- 
action |A|, ii) this does not imply pairing at all concen- 
trations, hi) we can predict if there is pairing or repulsion 
in a particular case by our theory. The cluster approach, 
however, has several limitations, the main size effect be- 
ing that W=0 pairs are possible at discrete values of the 
hole concentration. In Cu^Oi with 4 holes, we are doping 
with one electron, but in other cases we are far from the 
physical concentrations. However, we are not yet trying 
to make quantitative predictions, rather our point here 
is that of testing our approach against exact solutions, 
which is only feasible in small clusters. 

VII. FLUX QUANTIZATION AND PAIR 
SYMMETRY 

If a magnetic field is confined to a hole in any material 
(flux tube) the flux <j) is quantised in integer multiples of 
the fundamental quantum </>q — — ; a flux <j) = </>o can 
be gauged away, and any physical property, for example 
the ground state energy, is a periodic functions of <f> with 
period <po. 

Bulk superconductors quantize the flux through a 
hole in integer and half- integer multiples of 4>q, because 
the quasiparticles that screen the vector potential carry 
charge 2e. In finite systems the signature of supercon- 
ductivity is a ground state energy minimum at 4> — 
that is separated by a barrier from a second minimum at 
<f> = 4>o/2. With increasing the size of the system, the 
energy (or free energy, at finite temperature) barrier sep- 
arating the two minima becomes macroscopic, and bulk 
superconductors can swallow up only an integer or half 
integer number of flux quanta. As emphasized by Can- 
right and Girvin |27|] , the flux dependence of the ground 
state energy is definitely a most compelling way of test- 
ing for superconductivity, and the existence of the two 
minima separated by a barrier is a strong indication of 
superconducting flux quantisation. 

In Rcf. J27| , superconducting pairing was obtained by 
assuming a negative U; a ribbon shaped cluster was closed 
on itself with periodic boundary conditions along its 
length, and the flux was inserted in the hole. In the 
present problem, with a repulsive Hubbard model, the 
mechanism of attraction is driven by the C^ v symmetry, 
and cannot operate with such an unsymmetric geometry. 
The flux must be inserted in such a way that the system 
is not distorted. On the other hand, we cannot make 
holes in our small cluster because it would fall apart in 
disconnected pieces. One should consider larger clusters 
like CM13O36, which allow W=0 solutions for nh > 10, 
however the number of configurations > 10 12 is outside 
the scope of exact diagonalisation methods. 

So, we keep the CU5O4 cluster geometry, but modify 
its topology by adding a small hopping td between the 
external Cu's, in order to introduce a closed path around 
the centre, where screening currents can respond. Each 



td bond forms a closed triangular loop with the central 
Cu at the vertex (see Figure I). This geometry is a com- 
promise, because the magnetic field penetrates our small 
cluster; however, it lends itself to an extension to the 
full plane, such that only the 4 central plaquettes feel 
a magnetic field, and the rest of the plane only experi- 
ences a vector potential (see below) . Finally, we observe 
that a flux of the order of a fluxon in a macroscopic sys- 
tem would be a small perturbation; in the small cluster, 
however, the perturbation is small only if the hopping 
integral td is taken small compared to t. Numerically, 
the computations were performed with td = ±0.01 eV. 

We introduce a tube carrying flux <f> inside each of the 
triangles formed in this way. Every bond collects the 

2tH f A dr 

Peierls phase ; by symmetry, t is unaffected by 

the flux, while 

t d -» t d e^T (57) 

for a clockwise path, and the complex conjugate expres- 
sion a counterclockwise path. 

A. Superconducting flux quantisation: numerical 
results 

According to Table II, A^(n^) is negative and pairing 
results at <f> = for 7^=4,8 and 16; in all three cases 
we found that the ground state energy Eh{rih,4>) as a 
function of <j) has clearly separated minima at zero and 
half a flux quantum. Moreover, our criterion for pairing 
(A < 0) also leads us to a much more stringent crite- 
rion for superconducting flux quantisation than is drawn 
from the literature, since we need that both minima in 
the ground state versus flux curves also correspond to 
negative A. This is a much clearer signature of super- 
conducting flux quantisation than the generally accepted 
presence of the two minima, because it implies that the 
superconductor remains a superconductor after swallow- 
ing up the half flux quantum. Therefore, we computed 
Ah{nh,<p) m order to determine the flux dependence of 
the effective interaction. When A/j(n^,0) < 0, then 
Aft(n^, 4p) is also negative. For small enough t d , the 
response function 

R = &h(n h ,<j>) - A h (n h ,0) 

M 

is an intrinsic property of the original cluster with td — 
0. In Figure II we show R for several Uh values versus 
All the R curves have a local minimum at d> = 0, 
where they vanish; <f> = 0o/4 is a maximum and a second 
minimum occurs at <f> = 4>o/2; the = 4 curve is reduced 
by a factor of 3. The barrier gets lower with increasing 
nh, but the same qualitative trend can be seen in all 
cases. The numerical data also show that changing the 
sign of td produces a rigid shift of the curves by ^ 
such that the two minima interchange their places. 
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In Figure II we also report the absolute value 
| (*o(0) |*o(0)) I of the overlap between the the ground 
state eigenvectors in the presence and in the absence of 

the flux, for n h = 4 . It is clear that | ^ (%-)l*o(0)) | = 

0, and therefore the pairing state at zero flux and half 
fluxon are orthogonal. There is a clear analogy with the 
BCS theory; in that case, the Cooper wavefunction has 
s symmetry and the total magnetic quantum number of 
the pair vanishes in the absence of flux, but not at half 
a flux quantum p§| ]. Similar results for the overlap are 
obtained for the other nh values which correspond to par- 
tially filled shells. 

Our code automatically classifies the eigenvectors ac- 
cording to the IRREPS of C^ v . For positive td the point 
symmetry of the ground state wavefunction changes from 
X B 2 (x 2 - y 2 ) at 0=0 to 1 k 1 {x 1 + y 2 ) at <j) = For 
negative td the symmetry labels of the two minima are 
interchanged. For electron pairing, the symmetry of the 
states is the same as in the hole case. 

Since the vector potential lowers the symmetry, the 
eigenvectors cannot generally be classified according to 
the IRREPS of C± v \ however numerical data show that 
at half fluxon, the symmetry is dynamically enhanced 
(see below). 

B. Group theory aspects of superconducting flux 
quantisation 

These findings are required by general symmetry prin- 
ciples. In the absence of td, the full invariance Group 
of the cluster is S4 and the interacting ground state is 
degenerate, with ^-Ai and 1 I$2 components. A nonzero td 
at <fi — reduces the symmetry to the C± v subgroup; it 
turns out that with a positive td the expectation value of 
the magnetic perturbation is negative on 1 B 2 and posi- 
tive on 1 A\; therefore the ground state is 1 B 2 at td > 
but changes symmetry if the sign of td is reversed. Upon 
switching the vector potential A, the Cu-Cu hopping is 
complex and chiral, so the symmetry is lowered again 
from C4,, to its subgroup Z4, which contains only the ro- 
tations. Since is abelian, there are no degeneracies 
for a generic 4>, so there are no W=0 pairs and repulsion 
prevails. With increasing the flux from 0, the ground 
state energy increases to a maximum. Then it decreases 
because, at (j> = the Cu-Cu hopping of Equation j57| ) 
becomes —td, which is real; then the full C^ v symmetry 
is restored, resurrecting the W=0 pairs. The recovery of 
Ci v at <f> = % enables us to assign the eigenvectors to 
the IRREPs, as noted above. The change of symmetry 
of the pair is also readily understood: the perturbation 
caused by td > at <fi = becomes the opposite at half 
fluxon, so the 1 Ai state is lowest now. The signature 
of superconducting pairing is not only the existence of a 
well defined second minimum at half flux quantum, but 
also the fact that it corresponds to a A <0 situation, like 
at <j> = 0. 



This symmetry argument extends to the full plane. To 
see that, consider the pattern of Figure III. Here, the Cu 
sites are marked by X and the Oxygen sites by O; the 
black dots stand for tubes carrying flux </> each, symmet- 
rically disposed around the central Cu. Varying </> by an 
integer multiple of 0o corresponds to a gauge transfor- 
mation leaving all the physical properties invariant. The 
arrows help to visualise a convenient choice of the gauge 
at general <f>. Namely, running along an oriented bond in 
the sense of the arrow, 

J A-dr=^; (59) 

along the other Cu-0 bonds, not marked in the Figure, 
J A • dr = 0. One sees that in this way the flux through 
any closed path corresponds to the number of tubes sur- 
rounded by the path. The reflection operations of C± v 
are equivalent to </> — * — 0, reverse the directions of the 
arrows and for a generic 4> the symmetry Group reduces 
to Z4. However, at <fi = 4r the reversal of the magnetic 
field in the tubes corresponds to a jump by </>o, and this 
is equivalent to a gauge transformation: this implies that 
the symmetry Group gets larger, the new symmetry op- 
erations being reflections supplemented by a gauge trans- 
formation. Indeed, it follows from Equation ( |59| ) that the 
hopping parameter becomes it along the arrows, while it 
remains equal to t along the unmarked bonds of Figure 
III. Any reflection operation simply changes the signs of 
all the hoppings along the marked bonds. Now consider 
the unitary transformation S which changes the signs of 
all the Cu orbitals along both diagonal, except the cen- 
tral Cu. Since S also has the effect of reversing all the 
arrows, a x S is a symmetry, for all reflections a in C^ v . 
Moreover, since the product of two reflections is a rota- 
tion, the Group C& v including the rotations and the re- 
flections multiplied by S is isomorphic to C^ v . The W=0 
pairs appropriate for half a flux quantum must involve 
two holes belonging to the degenerate IRREP of C^ v . 
In this way, at <p — 2 ^ e symmetry is restored, 
allowing again for pairing and negative A. The W=0 
quasiparticles have just the correct symmetry properties 
in the presence of the vector potential to provide super- 
conducting flux quantisation in macroscopic systems. 

VIII. CONCLUSIONS 

We have examined the properties of the W=0 pairs 
by performing numerical diagonalizations of the C115O4 
cluster for various fillings. Some of these fillings are not 
representative of the concentrations that have been re- 
alised in the cuprates, but our theory depends on sym- 
metry and the concentration range to which it applies is 
wider than that obtained experimentally. We have shown 
that the effective interaction between the two holes in the 
W=0 pair can be obtained by computing A by exact di- 
agonalisation or, alternatively, by an analytical, recursive 
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canonical transformation; we have detailed the latter ap- 
proach, and derived a weak coupling approximation that 
agrees with the numerical results for A and with a previ- 
ous diagrammatic analysis. Pairing occurs when A < 0. 
An approximate symmetry under charge conjugation ex- 
ists leading to electron pairing as well as hole pairing in 
the sense defined in the Introduction. The numerical data 
confirm that when the filling is such that W=0 pairs are 
involved, A is small in absolute value, while the other 
fillings lead to strong repulsion. In one case, the W=0 
pair leads to a small repulsion, showing that the exis- 
tence of pairing is not a general property independent on 
filling. In all cases, we found that pairing or its absence 
can be reliably predicted by studying the behavior of the 
system at weak coupling, which supports the approxi- 
mations that we performed in a study of the full plane 
in Ref( J2(|) at last in some concentration ranges. We 
stress, however, that the instability of the Fermi liquid 
against a pairing interaction does not grant superconduc- 
tivity, since there is a competition with other order pa- 
rameters. Further investigations are necessary to analyse 
this point very close to half filling, where the antiferro- 
magnetic order prevails and the behavior could require a 
strong coupling analysis. Moreover, we expect that the 
renormalisation of the dispersion relation cannot be ne- 



glected (19 



In the C115O4 cluster the exact diagonalisation results 
show that W=0 pairs quantise flux in the superconduct- 
ing way. The ground states in presence of zero and half 
fluxon have different symmetries, like in BCS supercon- 
ductors. The superconducting flux quantisation property 
is due to the fact that the symmetry Group appropriate 
at half flux quantum in isomorphic with C^ Vl and this 
is not limited to small clusters, but general. Flux quan- 
tisation and pairing fit well together, both being conse- 
quences of the same symmetry principle. 
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FIGURE CAPTIONS 

FIGURE I The Cu 5 Oa cluster with 4 flux tubes (black dots) 
carrying flux <f>. X stands for Cu. The dotted lines represent the 
bonds providing a closed path around the centre 



FIGURE II Solid line (right scale): | (# W>)|*o(0)> |. The other 
lines (left scale) show the dimensionless response function R of 
Equation (Ba) for 71^=4,8 and 16. Note that for n^=4, 
A h (4, - A h (4,0) « 0.3t dl but A h (4, ^) is negative. 



FIGURE III Pattern of the vector potential A due to 4 flux 
tubes (black dots) carrying flux <f>. X stands for Cu.The line 
integral of A along each bond parallel to the arrow is ^ . 



TABLES 


Symmetry 


g e(eV) 


a\ 


1 -1.643 


e ,bi 


3 -0.43 


ai 


1 0. 


e,bi 


3 3.93 


a\ 


1 5.143 


TABLE I. 


One-hole levels of the CU5O4 cluster, with their 


symmetry labels, degeneracies g, and energy eigenvalues e 


with t = 1.3 eV, e v = 3.5 eV and Ed = 0. 


n h 


A h (meV) 


4 


-15.7 


6 


1469.2 


8 -10.85 


12 


43.72 


14 


1109.2 


16 


-25.47 



TABLE II. Exact diagonalisation results for Ah(nh) 
(meV), using U p —6 eV and Ud=5.3 eV. For rih=4, 8 and 16 
pairing takes place, and at uh=12 the repulsion is drastically 
reduced. For 71^=6 and 14 the W=0 pairs are not involved 
and the normal repulsion develops. 
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